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ABSTRACT: Diffusion-limited annihilation, A + A —>■ 0, and coalescence, A + A —>■ A, 
may both be exactly analyzed in one dimension. While the concentrations of A particles in 
the two processes bear a simple relation, the inter-particle distribution functions (IPDF) 
exhibit remarkable differences. However, the IPDF is known exactly only for the coales- 
cence process. We obtain the IPDF for the annihilation process, based on the Glauber 
spin approach and assuming that the IPDF's of nearest-particle pairs are statistically in- 
dependent. This assumption is supported by computer simulations. Our analysis sheds 
further light on the relationship between the annihilation and the coalescence models. 
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1. Introduction 

One-species annihilation, A + A — > 0, and coalescence, A + A — > A, have been extensively 
studied as basic prototypes of diffusion-limited reactions [1-21]. In one dimension, these 
simple models yield themselves to exact analysis: exact results for coalescence have been 
obtained, for example, with the method of inter-particle distribution functions (IPDF) [12- 
18], and annihilation has been analyzed exactly through a Glauber spin formalism [19-21] 
(and by several other techniques). 

The distribution of distances between nearest particles, i.e., the IPDF, plays a fun- 
damental role in the Smoluchowski theory of diffusion-limited reactions and it has been 
studied extensively for various reaction models [21-30]. In the case of one-species coales- 
cence in one dimension the IPDF is obtained exactly, as a byproduct of the IPDF method 
[15]. For one-species annihilation in one dimension, several aspects of the kinetics are 
known exactly, including the time dependence of the concentration of particles, but not 
the IPDF. 

In this paper, we compute the IPDF for the annihilation process from the (exactly 
known) spin-spin correlation function of the Glauber spin formalism. Our derivation is 
based on the assumption that the IPDF's of different particle pairs are statistically indepen- 
dent. Although this may seem a strong assumption, comparison to computer simulations 
data suggests that it may be exact. At the very least, our approach yields a remarkably 
good analytic approximation for the IPDF of the one-species annihilation model. 

The rest of this paper is organized as follows. In the next section we briefly summarize 
the exact results known for the annihilation and the coalescence processes. We also discuss 
the similarities between the IPDF method and the Glauber spin formalism. In Section 3, 
we develop a relation between the spin-spin correlation function in the Glauber method 
and the IPDF of the annihilation process. In Section 4, we apply this technique to the 
computation of the IPDF for annihilation, and for annihilation with input, and compare 
to results from numerical simulations. In the discussion, in Section 5, we reexamine the 
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conjecture of statistical independence. We demonstrate the need for proof by showing that 
the conjecture is wrong in the case of the closely related model of coalescence. 



2. Coalescence and Annihilation 

Our models are defined on a one-dimensional lattice with lattice spacing Ax. Each lattice 
site can be either empty, or occupied with exactly one particle. Particles hop to the 
nearest site to their right or left at rate D/(Ax) 2 . On long length and time scales this 
yields normal diffusion with diffusion coefficient D. When a particle hops onto a site which 
is already occupied, a reaction takes place: in the case of coalescence the impinging particle 
is removed, modeling the reaction A + A — > A, and in the case of annihilation both particles 
(the hopping particle and the target) are removed, modeling A + A — > 0. 

Coalescence has been analyzed exactly by the IPDF method (also known as the 
method of empty intervals) [see ref. 15 for a review]. In this approach, one defines E n (t) as 
the probability that an arbitrary sequence of n consecutive sites is empty at time t. From 
the E n one can compute various quantities of interest. For example, the concentration of 
particles is given by 

The E n satisfy a simple, closed master equation which can be solved exactly not only 
for the basic coalescence model, but also for a number of variations, including the back 
reaction, A — > A + A, input of A particles, inhomogeneous initial distributions, and finite- 
size lattices [11-18]. 

Although the analysis can be carried out in discrete form, it is much simpler (and 
more transparent) to consider the continuum limit. This is done by defining x = nAx and 
then letting Ax — > 0. The empty interval probabilities are replaced by the two-variable 
functions, E(x,t), and the concentration, for example, becomes 

dE(x,t) 
dx — 



The IPDF p(x,t), i.e., the probability that the nearest particle to an arbitrary particle is 
a distance x away at time t, can also be obtained from E(x, t): 

c(t) P (x,t) = ^^. (3) 

Throughout the remainder of this paper we will confine ourselves to the continuum limit. 

In the case of simple coalescence the equation of motion for the empty interval prob- 
ability is 

^ - 2D^ (4) 
dt ~ dx 2 ' [ } 

with boundary conditions 

E(0,t) = l and E(oo,t) = 0. (5) 

For generic homogeneous initial conditions (excluding exotic fractal distributions of the 
particles [11]) the system arrives at a universal long time asymptotic regime, where 



and 



x x 2 

P<yX '^ ADt eXp(y ~8Dt^ as t^oo. (7) 



The IPDF can be put in scaling form in terms of the dimensionless interparticle distance 
£ = c(t)x: 

7T 7T ^ 2 

p(£,t) = c(t)p(x,t) -> -^exp(--y), as ^ oo. (8) 

Notice that in the long time asymptotic limit it becomes stationary. 

The annihilation model is most easily analyzed in terms of its dual Glauber spin 
process (in the zero temperature limit) [19,21]. In the latter, each lattice site is occupied 
by a 'spin' variable which assumes one of two states, a = ±1. Any spin adjacent to a spin 
of opposite state may change its state, at rate D/(Ax) 2 . If we associate a particle with 
each pair of anti-parallel spins (i.e., the particles are at the boundaries between alternating 
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domains of parallel spins, Fig. 1), these particles are seen to diffuse and annihilate upon 
encounter, exactly as in the annihilation model. 

Consider the spin-spin correlation function, G(x' , x" ,t) = (cr(x', t)a(x", £)), where 
the angular brackets denote an average over different realizations of the process. For 
our purpose it is sufficient to limit the discussion to homogeneous initial distributions of 
'particles', in which case G(x' ', x" , t) is a function of x = x" — x' . The correlation function, 
G(x,t), satisfies the diffusion equation [19,21] 

dG d 2 G 

with boundary conditions 

G(0,t) = l and G(oo,t) = 0. (10) 

The concentration of particles is obtained from the correlation function, as 

l dG(x,t) 

«(*) = -2-^-U=o- (11) 

Notice that G satisfies the same equation as E, and for generic initial conditions it attains 
the same asymptotic limit as E. It follows that the concentration of particles for annihi- 
lation is exactly one half of the concentration in the coalescence process (in the long time 
asymptotic limit). In fact, the two processes maintain a similar simple relation at all times 
[3]. 

Another approach to the annihilation model relies on P e (x,t), the probability that 
there is an even number of particles in an arbitrary interval of length x at time t [20] . P e 
bears a simple relation to the spin-spin correlation function: G = 2P e — 1, and hence it 
too satisfies Eq. (9), but with P e (oo,t) = 1/2. The advantage of this approach is in that 
it deals with the particles directly and it does not require the dual process with spins, but 
it is essentially the same. 

In spite of these exact analyses, there is no obvious way to determine the IPDF 
for the annihilation process. Computer simulations reveal that the long distance tail of 
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the IPDF falls off exponentially, as e~ x [15,21]. This is also supported by non-rigorous 
theoretical arguments [15]. In this respect the annihilation process is radically different 
from coalescence (compare Eq. 8). In the following, we develop a method for calculating 
the IPDF in the annihilation process. 

3. Relation between IPDF and Correlation Functions 

For the annihilation process, the IPDF may be related to the spin-spin correlation function 
by assuming that at any particular time the particles (or the spin flips) constitute a renewal 
process on the line [31]. The IPDF is the the renewal probability density. For simplicity, 
we shall omit the time variable, since it plays no role in the spatial renewal process of the 
particles. 

Let p\\ (x) be the probability that two spins separated by a distance x are parallel. 
Then, the spin-spin correlation function is simply 

G(x)=p||-(l-p||) = 2p||(a;)-l. (12) 

Two spins would be parallel if there are no spin flips in between them, or if there is an even 
number of spin flips in between them. (Indeed, p\\ is identical to P e .) Let po(x) be the 
probability that the first flip from a given spin occurs at a distance x from it (the forward 
waiting distribution function), po is related to the renewal probability distribution (the 
IPDF, p) through [32]: 

f cx'p(x') dx> f°° , 
po(x) = / = c p(x ) dx , (13) 

J X X J X 

where c = 1/ J °° xp(x) dx is a normalization constant which happens to be equal to the 
concentration of particles. The probability that there are no spin flips up to a distance x 
from a given spin is 

p OO 

Po(s) = / Po(x')dx'. (14) 

J X 
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(Notice that p' is identical to the empty interval probability E, defined for the coalescence 
process.) Finally, let pi(x) be the probability that there are no spin flips up to a distance 
x from a given spin flip: 

Pl (x)= r P (x')dx' = ^. (15) 

Jx c 

We can now express p\\ as the sum of p' , and po convoluted with p an odd number 
of times (this yields an even number of spin flips) and convoluted with p\ . Working with 
the Laplace transform of the distribution functions (which we shall denote by a tilde), the 
convolutions assume the simpler form of products: 

P\\ = Po +Po(P + P 3 H )pi = 2TT^- ( 16 ) 

11 s s z 1 + p 

Here s is the Laplace transform variable (conjugate to x), and the last equality has been 
obtained using Eqs. (13)-(15) and standard properties of the Laplace transform. From 
Eq. (12), we finally obtain the relation between the spin-spin correlation function and the 
IPDF: 

*) = TW >.M^(i-sCM). (17) 

Assuming statistical independence of consecutive IPDF's one can also derive an ex- 
pression for the density-density correlation function c^ 2 \x,t), i.e., the joint probability of 
finding any two particles separated by a distance x at time t. Because there can be any 
integer number of particles between two arbitrary particles, the density-density correla- 
tion function is given by a sum over convolutions of the IPDF. In terms of the Laplace 
transform: 

=p + p 2 +f + --- = -^ = (18) 
c 1— p 2h(s) 

4. IPDF's for the Annihilation Process 

Let us now apply the above formalism to annihilation in the long time asymptotic limit. 
In this case, for generic initial conditions the correlation function approaches a universal 



form: 



x 

G(x,t) = erfc( 2 ^ / __ ) as ^ oo, (19) 



and the concentration reaches the asymptotic limit 



c(t) = OT' 88 ^ °°' (20) 

With these expressions, we get from Eq. (17), 

1 — ^pKkse k s erfc(fcs) 
^ 1 + v/yf/cs e fc2s2 erfc(/cs) ' 

where k = V2Dt. This can be inverted either numerically or through series expansions. 
For small x we find (in scaling form, with £ = cx) 

\p = <-\* 2 e + ^e + --- (22) 

The linear term has been derived exactly by Amar and Family [21]. For large x we get the 
asymptotic form (by contour integration) 

< 23 > 

where a = —zq ~ 0.357835 is the absolute value of the smallest negative root of 1 + 
\/7r^exp(2 2 )erfc(2) = 0. Thus, we do confirm the large distance exponential decay of the 
IPDF. 

We have performed a similar analysis for the transient regime at early times. Here, a 
discrete formulation is more appropriate, but that is derived in a straightforward, analogous 
way to the continuum limit of Section 3 [31]. Some typical results are shown in Figs. 2(a) 
and (b). In Fig. 2(c) we compare simulation results at a very late time, well after the 
long time asymptotic limit sets in, to the analytical IPDF obtained from Eq. (21). There 
is perfect agreement between theory and simulations, to the level of the statistical noise. 
Notice, however, that it is difficult to obtain good statistics in the long time asymptotic 
limit, due to the small number of surviving particles. 



We next analyze annihilation with a steady homogeneous input of particles at rate R 
per unit time per unit length. This model has the advantage of a stationary state, making 
it possible to obtain exceptionally good statistics for the IPDF from simulations. In this 
case, the spin-spin correlation function satisfies the equation [19] 



with the stationary solution 



(r = R/D). The stationary concentration is 

C " 2Ai(0) ' [ ' 

where the prime denotes differentiation of the Airy function. 

In Fig. 3 we compare simulation results to the analytic IPDF obtained from Eqs. (17), 
(25) and (26). Again, the agreement is excellent. The slight discrepancy at the peak of the 
curve may be a result of discreteness. Indeed, this discrepancy becomes smaller for smaller 
steady state concentrations. For short distances, we get (l/c)p(£) = 4(Ai(0)/|Ai'(0)|) 3 £. 
For long distances, we find an exponential decay of the IPDF tail, similar to pure annihila- 
tion. This is interesting, in view of the fact that the spin-spin correlation function decays 
as exp(— x 2 ) for pure annihilation, but as exp(— a; 3 / 2 ) in the case of input. 

5. Discussion 

Our analytical approach would be exact if the IPDF's for different particle pairs in the 
annihilation process were statistically independent. Only then could the convolutions 
leading to Eqs. (17) and (18) be justified. Therefore, the excellent agreement between 
simulations and the analytical computations leads us to the following conjecture: 

The distribution functions of the inter-particle distances of different particle pairs in the 
annihilation process are statistically independent from each other. 
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We assume, of course, that at time t = the initial distribution of particles is such 
that IPDF's are statistically independent. If such is the case, our results suggest that the 
statistical independence would be maintained at all times. 

It is not obvious why this conjecture should be true. In fact, we now show that in the 
coalescence process the IPDF's of different particle pairs are statistically dependent. The 
density-density correlation function for the coalescence process has been obtained exactly 
[33, 34]: 

c (2) (0 = c 2 [l -e-^ 2 + VK 2 erfc(^)]- (27) 

On the other hand, we can compute assuming statistical independence, from Eqs. (8) 
and (18). In spite of a rough qualitative agreement, there are important discrepancies: 
for example, the exact result of Eq. (27) shows that approaches the long distance 
limit as c 2 — ~ exp( — §£ 2 ), monotonously, while statistical independence predicts an 
exponential approach with oscillations, c 2 — ~ exp(— a£) cos(5£ + </>), (a, 6, and (p are 
known constants [34]). 

This poses the following puzzle. Consider diffusion-limited polymerization of n-mers, 
A n . The polymers diffuse on the line, with a diffusion constant independent of their size, 
and polymerize upon encounter: Ai + Aj — > A i+ j. Initially, all particles are monomers, 
A\. This process codes for coalescence and for annihilation simultaneously [10]: If we 
disregard size, and consider all polymers as equivalent particles, the process is analogous 
to coalescence, A + A — > A. On the other hand, if we focus only on polymers of an odd 
number of monomers, the process is equivalent to annihilation, A + A — > 0. Imagine now 
that we start this "master-process" and we look at a snapshot of the system at some later 
time. There are polymers of odd and even size, roughly in equal amounts, and apparently 
well mixed (there is no segregation into clusters of odd and even sized polymers). From our 
foregoing discussion we conclude that the distances between any pair of nearest polymers 
are statistically dependent (we are looking at coalescence). But if our conjecture is true, 
the distances between nearest polymers of odd size (annihilation) — a finite subset of all 
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polymers — are statistically independent! There is of course no contradiction, however, the 
spatial distribution of the polymers in the master- process must then be rather peculiar. 

In conclusion, we have presented strong numerical evidence that the distribution func- 
tions of distances between nearest particle pairs in the annihilation process, and in the 
annihilation process with input, are statistically independent. If we accept this as true, 
then the IPDF for the annihilation model has been computed exactly. Our computations 
reproduce the exact short distance limit of Amar and Family [21]. We also confirm the 
exponential decay of the IPDF tail. The IPDF for the case of annihilation with input is 
similar to that of pure annihilation, but the depletion zone near the origin is narrower. 
This is the result of the input, which tends to bring the IPDF closer to a pure exponential, 
because of its random nature. 

Our work suggests a most significant difference between annihilation and coalescence: 
the IPDF's are statistically independent for the former, but not for the latter. An interest- 
ing result is that the density-density correlation function exhibits oscillations in the case of 
annihilation, but not for coalescence. Such oscillations are typical of fluids with hard core 
repulsion interactions, and are expected because of the effective repulsion between particles 
due to the reactions. Their absence in the case of coalescence is a striking consequence of 
the correlations between IPDF's. 

The excellent agreement between our analytical derivations and simulations suggests 
that our approach is exact. We have not yet been able to prove the necessary conjecture 
of statistical independence (for the annihilation model). This remains an open problem. 
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CAPTIONS 



Figure 1: Relation between the Glauber spin model and the annihilation process. The 
A particles (bottom) correspond to the spin flips, or domain boundaries, in the spin 
process (top). 

Figure 2: Comparison of simulated IPDF's to analytical derivations for the annihilation 
process, A + A — > 0. Simulation results are shown in histogram form for (a) t=l, 
(b) t = 10, and (c) t = 1000. All simulations began with a homogeneous density of 
particles, 1/2. 

Figure 3: Comparison of the simulated IPDF (circles) to the analytical derivation (solid 
curve) for the annihilation process with input. The input rate (R = 0.0001) is low 
enough so that the continuum limit applies (c ss 0.0213 <C 1). Simulation results 
represent an average over 16 x 10 6 time units (performed after the steady state was 
reached). 
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